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We investigate the stability of cosmological scaling solutions describing a barotropic fluid with 
' p = (7 — l)p and a non- interacting scalar field (f> with an exponential potential V(4>) = Voe~ K<1 '. 

, We study homogeneous and isotropic spacetimes with non-zero spatial curvature and find three 

■ possible asymptotic future attractors in an ever-expanding universe. One is the zero-curvature 

power-law inflation solution where fi^ = 1 (7 < 2/3, ft 2 < 37 and 7 > 2/3, k 2 < 2). Another is 
the zero-curvature scaling solution, first identified by Wetterich, where the energy density of the 
scalar field is proportional to that of matter with 0,$ = S'j/k 2 (7 < 2/3, k 2 > 37). We find that this 
matter scaling solution is unstable to curvature perturbations for 7 > 2/3. The third possible future 
asymptotic attractor is a solution with negative spatial curvature where the scalar field energy 
density remains proportional to the curvature with Q.^ = 2/k 2 (7 > 2/3, k 2 > 2). We find that 
solutions with = are never late-time attractors. 
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O I. INTRODUCTION 



There have been a number of studies of spatially homogeneous scalar field cosmological models with an exponential 
potential, with particular emphasis on the possible existence of inflation in such models if the potential is sufficiently 
flat . Although modern supergravity theories generally predict exponential potentials that are too steep to drive 
■ inflation, these models still have important cosmological consequences. For example, in models with barotropic matter, 
rT\ such as dust or radiation, there exist spatially homogeneous and isotropic 'scaling solutions' in which the scalar field 
energy density tracks that of the perfect fluid |6|-^) , so that a significant fraction of the energy density of the Universe 
at late times may be contained in the homogeneous scalar field whose dynamical effects mimic the barotropic matter. 
The tightest constraint on these cosmological models comes from primordial nucleosynthesis bounds on any such relic 
density during the radiation dominated era ||. More recently attention has been directed to the possible effect of 
such a scalar field on the growth of large-scale structure in the universe [p|-|r2|] . 

A phase-plane analysis of the spatially homogeneous and isotropic zero curvature models S has shown that these 
scaling solutions are the unique late-time attractors whenever they exist. The stability of these scaling solutions in 
more general spatially homogeneous cosmological models was studied in |l3). In this article we shall study scaling 
solutions in the general class of spatially homogeneous and isotropic models with non-zero curvature. 
The governing equations for a self-interacting scalar field with an exponential potential energy density 

V = Voe~ K t (1.1) 

where Vq and K are positive constants, evolving in a Robertson- Walker spacetime containing a separately conserved 
perfect fluid, are given by 

H = -\hp- l + 4> 2 )-K, (1.2) 
Pi = -3jHp 1} (1.3) 
4> = -3H^+kV, (1.4) 

subject to the Friedmann constraint 
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H 2 = l -( Pl + l -tf + V)+K, 



(1.5) 



where K = —kR~ 2 and k is a constant that can be scaled to 0, ±1, H = R/R is the Hubble parameter, an overdot 
denotes ordinary differentiation with respect to time t, and units have been chosen so that 8irG = 1. In the above we 
have assumed that the perfect fluid satisfies the barotropic equation of state 



P-y = (7 - !)P7' 



(1.6) 



where 7 is a constant which satisfies < 7 < 2. We also note that the total energy density of the scalar field is given 
by 
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V. 



Defining 



V 



and the new logarithmic time variable r by 



V3H 



dT -H 



3iJ 2 ' 



equations (1.2-1.4) can be written as the three-dimensional autonomous system: 



/3 2 3 
x = -3x + W -ny + -x 



7 - - )n+-(l + 2x z -y z ) 



y =- 2 v 



n' = an 



-nx 



7-- )n+-(l + 2x z -y J 



7- 



(n-l) + -(2x 2 -y 2 ) 



where a prime denotes differentiation with respect to r, and equation (1.5) becomes 

1-fl-x 2 -y 2 = KH~ 2 , 

where 



n 



3H 2 



2 , 2 
x +y . 



(1.7) 

(1.8) 

(1.9) 

(1.10) 
(1.11) 
(1.12) 

(1.13) 
(1.14) 



II. QUALITATIVE ANALYSIS 
A. Invariant Sets and Monotone Functions 



The physical region of the state-space is co nstrained by the requirement that the energy density be non-negative; 
i.e., O > 0. Furthermore, from equation ( 1.13 ) we find that in the variables used the state-space is bounded for k = 
and k = — 1 but not for k = +1. 

Geometrically the zero-curvature models (k = 0) are represented by the paraboloid 

n + x 2 + y 2 = 1, (2.1) 

in the (f2, x, y) state-space. It is therefore possible to divide the complete state-space up into a number of invariant 
sets. (See Table |.) 



2 



The existence of a monotone function in any invariant set rules out the existence of periodic orbits, recurrent 
orbits, and equilibrium points in that set and serves to determine the asymptotic behaviour in that set jl4 15 1. In the 
three-dimensional sets A and C we define the function 



Z = 



n 2 



From Eqns. (1.10-1.12) we find that 



(n + x 2 + y 2 -l) 2 ' 
Z' = 2(2 - 3j)Z. 



(2.2) 



(2.3) 



For 7 > 2/3 we have that Z is a monotone decreasing function along orbits in sets A and C, which therefore implies 
that r2 — > along these orbits. For 7 < 2/3 we have that Z is a monotone increasing function along orbits in sets A 
and C which also implies that K/QH 2 — ► along these orbits. (The case 7 = 2/3 will be dealt with separately in Sec. 
ETC.) Hence we observe that in the non-zero-curvature models either the curvature approaches zero or the energy 



density of the fluid approaches zero as time evolves forward. That is, the asymptotic behaviour of these models can 
be completely determined by analyzing the zero-curvature models (see Q), and the scalar field models with no fluid 
matter (see 



B. Local Stability of the Equilibrium Points 



Due to the existence of the monotonic function Eq. (2.2), all we need to do is locally analyze the equilibrium points 
of the system (1.10-1.12). There are seven equilibrium points, labelled P1-P7. The equilibrium points and their local 
stability are listed in Table |l[ 

The exact solution corresponding to the equilibrium point PI is the standard Milne model. Equilibrium points 
P2 and P3 correspond to a massless scalar field models which are essentially equivalent to stiff perfect fluid 
Friedmann-Robertson- Walker models. The remaining equilibrium points and their corresponding exact solutions 
are discussed later in the text. 



C. The Bifurcation Values 



Bifurcations of the the system ( |1.10| - |l.l2| ) occur for various values of the equation of state parameter for the perfect 
fluid, 7, and the scalar field potential parameter, n. The bifurcation values are 7 = 2/3, k 2 = 2, k 2 = 6, and k 2 = 3j. 
Each of these bifurcation values will be discussed in turn. 



1. 7 = 2/3 



If 7 = 2/3 then the equilibrium points of the system ( 1.1C - L12] ) are points P2, P3, P4 as well as the non-isolated 
lines of equilibria given by 

LI = (fi = fio, x = 0,y = 0) 



and 



L2 = (n = n a 



V2_ 

\fin V3k 



where < flo, ft s < 00. These two lines of equilibria are the degeneracies of equilibrium points PI, P5, P6 and P7. 
Fortunat ely the analysis of the dynamics of the system are simplified through the observati on that th e function defined 
by Eq. ( |2.2| ) is constant when 7 = 2/3. This implies that the dynamics of the system (|1.10| - |lT2|) are restricted to 
level surfaces of the function Z = Z(x, y, £1). Essentially this implies that dynamics of the three-dimensional system 



( l. 10 - 1.12 ) can be thought of as a one-parameter family of two-dimensional surfaces (paraboloids in this case), where 
the dynamics on each surface are identical and equivalent to the dynamics on any other surface, (including the surfaces 
represented by Z = and l/Z = 0). The future asymptotic attractor for all models is the power-law inflationary 
model represented by equilibrium point P4 if k 2 < 2. If k 2 > 2 then the matter scaling solutions represented by the 
line of equilibria L2 are the future asymptotic attractors for all models. 
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2. k 2 = 2 



If k 2 = 2 then points P4 and P5 coalesce. As k 2 — * 2~ , P5 ^P4, and as k 2 increases past 2, the stability of 
one of the eigendirections changes. In essence P4 and P5 undergo a transcritical bifurcation ]l6|j . The point P4 is 
a saddle node when k 2 = 2. If k 2 = 2 and 7 > 2/3 then both the negative and zero curvature models tend to the 
equilibrium point P4. For the positively curved models, this same point acts like a saddle with a one-dimensional 
unstable manifold when k 2 = 2 and 7>2/3. If k 2 = 2 and 7 < 2/3 then all models are attracted to the matter 
scaling solution represented by P7. 



If k 2 = 6, then equilibrium points P4 and P2 coalesce. However, the future asymptotic behaviour of the system 
can be determined directly from Table |l[ 

4. K 2 = 37 

If k 2 — 37 then P7 and P4 coalesce. If 7 < 2/3 then this point attracts those orbits in the physical state space, i.e., 
it is the future attractor. If 7 > 2/3 then there exists a one-dimensional unstable manifold. If k 2 = 3j and 7 > 2/3 
then the negatively curved models are attracted to the point P5, the zero curvature models are attracted to the point 
P4. 

D. Zero-Curvature Models 

The zero-curvature models are contained in the two-dimensional invariant set B U £. A qualitative analysis of 



this plane-autonomous system is given in 
corresponds to the equilibrium point P4, w 



The well-known power-law inflationary solution for k < 2 |2j-|fj] 

mich is shown to be stable (i.e., attracting in the two-dimensional invariant 
set BU £) for k 2 < 37 in the presence of a barotropic fluid. 

In addition, for < 7 < 2 there exists a matter scaling solution corresponding to the equilibrium point P7, 
whenever k 2 > 37 The equilibrium point is stable in the two-dimensional invariant set B \J £, (a spiral for 
k 2 > 2Aj 2 /(9j — 2), otherwise a node) so that the corresponding cosmological solution is a late-time attractor in the 
class of zero-curvature models in which neither the scalar-field nor the perfect fluid dominates the evolution and we 
have 

n* = § • ( 2 - 4 ) 

The effective equation of state for the scalar field is given by 

= {P± + P<t>) 2x 2 
^ — 2 1 2 

which is the same as the equation of state parameter for the perfect fluid at this equilibrium point; i.e., 7^ = 7. The 
solution is referred to as a matter scaling solution since the energy density of the scalar field remains proportional 
to that of the barotropic perfect fluid according to py/pd, = (k 2 — 37)/37 ||. Since this matter scaling solution 
corresponds to an equilibrium point of the system (1.1C- 1T2] ), we note that it is a self-similar cosmological model 

& 

E. Non-Zero Curvature Models 

1. 7 > 2/3 

If k 2 < 2, then the future asymptotic state for the negative-curvature models and a subset of the positive-curvature 
models is the standard zero-curvature power-law inflationary solution, represented by equilibrium point P4, where 
Qs = 1. Previous analysis has shown that this power-law inflationary solution is a global attractor in spatially 
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homogeneous models in the absence of a perfect fluid (except for a subclass of Bianchi type IX models which recollapse) 
|0|. However, more recent analysis has shown that this power-law inflationary solution is also a global attractor in 
spatially homogeneous Bianchi class B models (as classified by Ellis and MacCallum |T^]) in the presence of a perfect 
fluid |18| with the same restrictions on 7 and k as above. 

If nr> 2, then the negative-curvature models are asymptotic towards the curvature scaling solution, represented 
by equilibrium point P5, where 

«0 = \ ■ (2.5) 

Part of the motivation here is to study the stability of these scaling solutions. It is known that the equilibrium point 
P5 is the future asymptote for the Bianchi type V and Bianchi type VII^ models (see 13 1§|]). 



2. 7 < 2/3 

If k 2 < 37, then power-law inflation, represented by equilibrium point P4, is again the future asymptotic state 
for the negative-curvature models and a subset of the positive-curvature models. But for k 2 > 37 the matter scaling 
solution, represented by equilib rium point P7, takes over as the future asymptotic attractor. 



The linearization of system ( 1.1C - L12] ) about the equilibrium point P7 yields two negative real eigenvalues and 



the eigenvalue (37 — 2). Hence the matter scaling solution is only stable for 7 < |. For 7 > | the equilibrium point 
P7 is a saddle with a two-dimensional stable manifold (lying in the set B U £ representing the zero-curvature models) 
and a one-dimensional unstable manifold. 

III. COSMOLOGICAL IMPLICATIONS 

In order to solve the flatness problem of the standard model, a period of accelerated expansion where the the 
spatial curvature is driven to zero (K/H 2 — » 0) is desirable jis|| . Power-law models of inflation driven by scalar fields 
with exponential potentials with k 2 < 2 provide an interesting model of inflation where exact analytic solutions are 
possible, not only for the background homogeneous fields jl]], but also for inhomogeneous linear perturbations [p0| . 
We have shown that the power-law inflation solution, represented by equilibrium point P4, with fi^ = 1 and zero 
curvature, is a stable late time attractor solution in models with spatial curvature and non-zero density of barotropic 
fluid for 7 > 2/3 and k 2 < 2 and for 7 < 2/3 and k 2 < 3j. 

Power-law inflation can also be driven by a barotropic fluid with 7 < 2/3. However, we have shown that this fluid 
dominated solution, corresponding to equilibrium point P6, is never stable for 7 > in the presence of a scalar field 
with an exponential potential. Instead, for 7 < 2/3 and k 2 > 37 we have shown that the matter scaling solution, 
corresponding to equilibrium point P7 with zero spatial curvature, is the late time attractor for all models with zero 
or negative curvature, and for a subset of models with positive curvature. 

This leaves a possible relic density problem Although the spatial curvature may become negligible, the energy 
density in a scalar field with an exponential potential cannot be diluted away by inflation driven by barotropic matter 
unless 7 — > 0. In typical slow-roll inflation models 7 ~ 0.1 plfl , and thus the exponential potential can have a 
non-negligible density during and after inflation. During the subsequent radiation dominated era (7 = 4/3) the scalar 
field rapidly approaches the matter scaling solution where f2^ = 4/k 2 . Models of primordial nucleosynthesis require 
fl^ < 0.13 to 0.2 (|). Taking the more conservative upper limit we can conclude that the existence of scalar fields with 
exponential potentials is incompatible with standard models of nucleosynthesis for k 2 < 20, unless the scalar field has 
for some reason not reached its scaling solution. Conventional inflation models with < 7 < 2/3 cannot prevent 
reaching its equilibrium value soon after inflation. 

Observational evidence suggests we live in a Universe with a present day energy-density Qq bounded between 0.1 
and 0.3 |22j |. An intriguing possibility is that the dominant dark matter in the Universe could be in the form of a 
scaling scalar field with an exponential potential, which would be compatible with a low density, but spatially flat 
universe for k 2 ~ 4. However, this would only be compatible with the nucleosynthesis bounds quoted above if the 
scalar field is far from its scaling solution at the time of nucleosynthesis 0. Such a model is also likely to be in 
conflict with conventional models of structure formation as an inhomogeneous scalar field behaves differently from 
the barotropic matter and scalar field gradients can exert large pressures that resist gravitational collapse on small 
scales [|l2]|l0 



We have shown that the matter scaling solution is always unstable to curvature perturbations in the presence of 
ordinary matter (7 > 1); i.e., the matter scaling solution is no longer a late-time attractor in this case. However, it 
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does still correspond to an equilibrium point in the governing autonomous system of ordinary differential equations 
and hence there are cosmological models that can spend an arbitrarily long time 'close' to this solution. Indeed the 
curvature of the Universe is presently constrained to be small by cosmological observations, so it is possible that the 
matter scaling solution could be important in the description of our actual Universe. Negative curvature spaces are 
interesting since they can support a variety of possible topologies, leading to, for example, hyperbolic null geodesic 
motions which can produce interesting patterns in the cosmic microwave background radiation p3| . At late times, all 
models with n 2 > 2, 7 > 2/3 and negative spatial curvature approach the curvature scaling solution, corresponding 
to equilibrium point P5. 

Classical cosmological tests such as number counts, or the luminosity of standard candles, should in principle be 
able to distinguish between the different possible low-density models. The deceleration parameter at any point in the 
phase-space is given by 

<b = -7P= (^-1)0 + 2^-^. (3.1) 

This reduces to the standard result go = 37/2 — 1 for the matter scaling solution, or go — for the curvature scaling 
solution. The dynamics of the cosmological scale factor in these these solutions is indistinguishable from the matter 
dominated or curvature dominated solutions. However the luminosity distance as a function of redshift depends not 
only upon the evolution of the scale factor, but also upon the spatial geometry, and hence the curvature scaling 
solution can in principle be distinguished from the curvature dominated Milne model. 

Recent results from high redshift supernovae searches suggest that the deceleration parameter is negative at present; 
i.e., the Universe is in fact accelerating pi} . Although it is possible to obtain go < when y 2 > 2x 2 , we have found 
from numerical experimentation that this is unlikely during a transition from matter scaling to curvature scaling 
regime when n 2 > 2. Instead it would only seem to be compatible with either barotropic matter with 7 < 2/3 or a 
scalar field dominated solution, corresponding to equilibrium point P4, with k 2 < 2. 



IV. CONCLUSIONS 



A complete qualitative analysis of the dynamical system describing the evolution of spatially homogeneous and 
isotropic models containing a non-interacting perfect fluid and a scalar field with an exponential potential has yielded 
the following results: 

1. The past asymptotic state for all expanding models is the massless scalar field solution. 

2. The future asymptotic state depends upon the values of the equation of state parameter for the fluid, 7, and 
the steepness of the potential n. The future asymptotic state is either 

(a) the standard power-law inflationary solution (if 7 < 2/3, k 2 < 37 or 7 > 2/3, k 2 < 2), 

(b) a matter scaling solution (if 7 < 2/3, k 2 > 37), or 

(c) in the case of the negative-curvature models, a curvature scaling model with no barotropic matter (7 > 2/3, 
k 2 >2). 

3. In non-zero curvature models with ordinary matter, (i.e., 7 > 1), the matter scaling solution is never stable. 

Current cosmological observations seem to indicate < 1 in ordinary matter which could be due to the presence 
of (negative) spatial curvature, and or some other form of energy density such as a scalar field with exponential 
potential. It is known that the matter scaling solutions are late-time attractors (i.e., stable) in the subclass of zero- 
curvature isotropic models S . It is also known that the matter scaling solutions are stable (to shear and curvature 
perturbations) in generic anisotropic Bianchi models when 7 < 2/3 J13| . However, when 7 > 2/3, and particularly 
for ordinary matter with 7 > 1, the matter scaling solutions are unstable; essentially they are unstable to curvature 
perturbations, although they are stable to shear perturbations Even though it is unstable, this matter scaling 
solution may still play an important role in describing the intermediate (or transient) physics between the past and 
future asymptotic attractors. 

The negative curvature scaling solutions may also be of importance since they are self-similar cosmological models 
corresponding to an equilibrium point of the dynamical system ( l.lOj - 1.12 ) where the scalar field still has a non- 



vanishing contribution to the energy density. Cosmological tests such as the magnitude-redshift relation for standard 
candles or detailed observations of the microwave background should in principle be able to determine whether these 
solutions could describe the present state of the universe. 
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TABLE I. Invariant sets for the dynamical system given by equations 



Label Set (Cl,x,y) Dimension Description 

34 £1 > and Q + x z + y 2 > 1 3 positive curvature, non-vacuum 

B Q > and Q + x 2 + y 2 — 1 2 zero curvature, non-vacuum 

C £1 > and Q + x 2 -f- y 2 < 1 3 negative curvature, non-vacuum 

T> Q — and Q + x 2 + y 2 > 1 2 positive curvature, no fluid matter 

£ O — and Q + x 2 + y 2 — 1 1 zero curvature, no fluid matter 

J 7 Q — and Q + x 2 + y 2 < 1 2 negative curvature, no fluid matter 



Liq i.i2). 
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TABLE II. Equilibrium points of the dynamical system (1. 10-1. 12) and values of the parameters 7 and k for which the 
equilibrium point is a future attractor. 



Label 


a 


X 


y 


3 R 


Eigenvalues 


Stability Conditions 


PI 











-1 


2-37 
-2 
1 


Never 


P2 





1 








6-37 
4 


Never 


P3 





-1 








6-37 
4 

3 + ^Sk 


Never 


P4 





K 
fa 

vo 


V D 





k 2 - 37 
k 2 - 2 


7 < §, k a < 37 
and 

7 > f , k 2 < 2 


P5 





V2 


2 


2 - k 2 
2 


2-37 


7 > §, k 2 > 2 


P6 


1 











37-2 
i (3 7 - 6) 


Never 


P7 


K 2 


2k 


^12 7 -6 7 2 
2k 





37-2 

f [(7 - 2) ± ^(7 - 2) 2 + 2J( 7 - 2) ( K 2 - 3 7 )J 


7 < §, k 2 > 37 
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